A heterogeneous zero-range process related to a two-dimensional walk model 
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We have considered a disordered driven-diflusive system defined on a ring. This system can be 
mapped onto a heterogeneous zero-range process. We have shown that the grand-canonical partition 
function of this process can be obtained using a matrix product formahsm and that it is exactly 
equal to the partition function of a two-dimensional walk model. The canonical partition function 
of this process is also calculated. Two simple examples are presented in order to confirm the results. 
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I. INTRODUCTION 



A Zero-Range Process (ZRP), which was first intro- 
duced in mathematical literature [1], is a widely studied 
Markov chain in statistical mechanics. A homogenous 
ZRP is a system of interacting particles which hop be- 
tween the lattice sites (or boxes) of a homogeneous lattice 
with the rates that depend only on the number of parti- 
cles in the departure site. This can be generalized to a 
Heterogeneous ZRP (HZRP) in which the particles hop 
between the lattice sites (or boxes) of a heterogeneous 
lattice so that the particles' hopping rates also depend 
on the departure site. Different types of these processes 
have been introduced and studied in related literature 
(for a review see and references therein). They have 
also a variety of applications in different fields of physics 
for instance in transport processes including traffic flows, 
shaken granular gases, particle condensation and phase 
separation It is long known and well-stablished 

that there is possibility for mapping a one-dimensional 
driven-difltusive system defined on a discrete lattice with 
periodic boundary conditions onto a ZRP [2I]. 
Recently the study of the one-dimensional driven- 
diffusive systems which can be mapped onto the walk 
models has drawn physicists attention Q-flll- It has 
been shown that the simplest system of this type, that is 
the Asymmetric Simple Exclusion Process (ASEP) with 
open boundaries, can be mapped onto a walk model con- 
sisting of different Dyck paths Q-dl- The ASEP is a sys- 
tem of hard-core particles moving in a preferred direction 
on a one-dimensional discrete lattice [12]. The hard-core 
particles are injected from the left boundary of the lattice 
and hop exclusively toward the right boundary. Finally 
they leave the lattice from the right boundary. It has 
also been shown that the Partially ASEP (PASEP) in 
which the hard-core particles can hop both to their left- 
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most and rightmost neighboring empty lattice sites can 
be mapped onto a weighted walk model f9|-[l0j. 
In a recent paper [ll| it has been shown that a driven- 
diffusive system consisting of two different classes of par- 
ticles defined on a discrete lattice with periodic boundary 
condition can be mapped onto a two-dimensional walk 
model. The system is defined on a ring of length N . It 
consists of M — 1 first-class particles and a single second- 
class particle sometimes called the impurity. Both the 
first and second-class particles hop in the same preferred 
direction, and there is no overtaking. The steady-state 
of this system was first studied using a matrix product 
method (for a review of the matrix product method see 
p^ ) in [l3|- It was also shown that the system can 
be mapped onto a HZRP if one considers the particles 
(empty lattice sites) in the driven-diffusive system as 
boxes (particles) in the HZRP [l^. It turns out that 
if the single second-class particle moves slow enough, a 
large number of empty lattice sites pile-up in front of it 
otherwise the number of empty lattice sites in front of 
the second-class particle is of order 1. In the equivalent 
HZRP these two phases correspond to the cases where 
the number of particles in the box associated with the 
second-class particle is either of order of the number of 
empty lattice sites or 1. In [ll| the author has shown 
that the grand-canonical partition function of this disor- 
dered driven-diffusive system is equal to that of a walk 
model defined on the first quarter of a two-dimensional 
plane. The random walker starts from the origin and 
takes consecutive steps either upward or downward, pro- 
vided that it does not lie on the horizontal axis. When- 
ever the random walker lies on the horizontal axis, it can 
either goes one step upward or moves horizontally. After 
taking iV — 1 steps it can get to any height between 
and iV — 1. The paths taken by the random walker are 
weighted. The grand-canonical partition function of walk 
model is equal to the sum of these weights. If one consid- 
ers the paths which contain exactly N — M upward steps, 
then the resulting partition function is exactly that of the 
disordered driven-diffusive system explained in [l^ | and 
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In this paper we start with a generaUzed disordered 
driven-diffusive system with long-range interactions de- 
fined on a one-dimensional discrete lattice with periodic 
boundary conditions. We show that the steady-state of 
this system, which is a product measure, can be obtained 
using a matrix product method. The matrices associated 
with the particles and the empty lattice sites satisfy an 
algebra which has not been studied thus far. We discuss 
how the driven-diffusive system can be mapped onto a 
HZRP. Both the canonical and grand-canonical partition 
functions of the HZRP are calculated. An equivalent 
walk model is also introduced and studied in detail. Fi- 
nally, two simple examples are provided to confirm the 
results. 



II. THE DISORDERED DRIVEN-DIFFUSIVE 
SYSTEM 

In the author has introduced a disorder driven- 
diffusive system consisting of M different particles which 
hop in a preferred direction on a closed ring of length N. 
The particles are labeled with /i where fj, — 1, - ■ ■ , M. A 
particle of type fj, hops to its rightmost neighboring lattice 
site with the rate p^, provided that the target lattice site 
is not already occupied, according to the following rule: 

fi — > /i with the rate p^. 

Note that the particles do not overtake; therefore, the 
sequence of particles is preserved. On the other hand, the 
particle's hopping rate is a constant and depends only on 
the particle type /i. Let us denote the number of empty 
lattice sites in front of a particle of type jjL as, n^. In the 
steady-state the probability of finding the system in the 
configuration {ni, ■ ■ ■ , riAf } is given by: 

P({ni,--- ,nM}) = -J—Tr{DiE''' ■■■DME^'^n (1) 

in which the normalization factor Z^si.m is called the 
canonical partition function of the system. One should 
also note that "■/i = N — M. In ([1]) the matrices 

and E stand for the presence of a particle of type /i and 
an empty lattice site respectively. It has been shown that 
these matrices satisfy the following quadratic algebra: 

D^E=—D^ , = !,■■■ ,M. (2) 

In the authors have introduced a matrix representa- 
tion for the quadratic algebra as follows: 

oo oo 

D,=Y,U{^)\0){^\ , = + (3) 

i=0 i=0 

in which f^{i) = n j; = Note that {\i)} is 

the standard basis of an infinite-dimensional vector space 



that is \i)j = Si J for i,j = 0,1,-- - ,oo. The steady- 
state probability distribution ([1]) can now be calculated 
and one finds: 

M 

P({ni,--- ,nM}) = ^ II fM (4) 

in which the normalization factor Z^^m is given by: 

M 

Zn,M ^ 2^ ^E^"^, (JV-M) 

Let us now consider a generalized disordered driven- 
diffusive system consisting of different types of hard-core 
particles. We assume that a particle of type hops one 
lattice site forward according to the following rule: 

fi p fi' — > /i p -J ■ /i' with the rate u^(rt^i). 

This is a disordered driven-diffusive system with long- 
range interactions since the hopping rate of the parti- 
cle of type /i, that is itp(n^t), depends on the number of 
empty lattice sites in front of it. The steady-state prob- 
ability distribution of this generalized disordered system 
P({ni, • ■ • ,nA/}) is a product measure and can still be 
written as a trace of a matrix product similar to ([T]). It 
is not surprising since this generalized disordered driven- 
diffusive system can be mapped onto a ZRP. Therefore 
it has a product measure steady-state. In order to verify 
this, note that the matrices D^j, and i? in ([T]) satisfy the 
following algebra: 

D^.E'^-D,, = f^{n,,)D^, ioT ti,fi' = !,■■■ ,M (6) 

provided that we define ffj,{i) — n}=i TTIJ) that 
/^(O) = 1. It is easy to verify that the matrix repre- 
sentation (jS)) satisfies the algebra (jG)). Using the algebra 
©J straightforward calculations show that the steady- 
state probability distribution of the generalized disor- 
dered driven-diffusive system ([T]) can be simplified to (|4]) . 
It should be noted that for the special case M^(n^) = 
(equivalently f^{n^) = — ^) the matrix representation 

([3]) satisfies both ^ and (H]); however, for a general 
u^{n^) the matrix representation (|3]) only satisfies ([6]). 

III. THE HETEROGENEOUS ZERO-RANGE 
PROCESS 

The generalized disordered driven-diffusive system in- 
troduced in the previous section can be mapped onto a 
HZRP. As for this mapping a particle of type ^ in the 
driven-diffusive system should be considered as the ^'th 
box in the ZRP. The number of empty lattice sites in 
front of that particle in the driven-diffusive system is 
then equal to the number of particles in the box asso- 
ciated with that particle. In this case if the number of 



3 




1 2 3 4 5 



FIG. 1: Mapping between a disordered driven-diffusive sys- 
tem (above) and a heterogeneous zero-range process (below) 
with a periodic boundary condition. 

particles in the /z'th box is then the rate at which a 
particle moves from /Lt'th box to {fi — l)'th box is given by 
M^(rip). Note that if the particles in the driven-diffusive 
system move forward, they move backward in the ZRP 
(see FIG. [ij . It can be seen the particle's hopping rate 
from one box to the next box depends on both the box 
number and the number of particles inside it, which can 
be at most of order of the total number of empty lattice 
sites in the driven-diffusive system. In FIG. [l]we have 
provided a simple sketch of such mapping. The steady- 
state of the HZRP defined above is still given by (|4]). 
This has already been obtained in [2] using a different 
ap pro ach. 

In [la] the authors have studied a homogeneous ZRP and 
calculated its steady-state probability distribution using 
a matrix product method similar to that which we have 
done here. Their algebra and its matrix representation 
are special cases of ^ and ([3]). The main difference 
is that in a homogeneous system the particles' hopping 
rates do not depend on the box number. 



IV. THE WALK MODEL 

In the author has shown that a simple disordered 
driven-diffusive system, consisting of a finite number of 
ordinary particles in the presence of a single impurity, 
which can be mapped onto a HZRP, has a walk model 
interpretation. In this section we aim to show that the 
generalized HZRP introduced in the previous section has 
a simple equivalent walk model in a special case. Note 
that if the particles' hopping rates are constant and do 
not depend on the number of empty lattice sites in front 
of each particle, the model converges into the process in- 
troduced and discussed in [l^l and p^ . 
Let us start with the definition of the disordered driven- 
diffusive system. We assume that the driven-diffusive 
system is defined on a ring of length N and consists of 
two different types of particles. There is one second-class 
particle (let us consider it as the first particle) which 
hops to its rightmost vacant neighboring lattice site with 
the rate U2{n) provided that there are n empty lattice 
sites in front of it. There are also M — 1 first-class par- 



ticles on the ring which are considered as 2nd, 3rd, ■ • • 
and Afth particles. Each first-class particle hops to its 
rightmost vacant neighboring lattice site with the rate 
ui{n) provided that there are n empty lattice sites in 
front of it. A simple example is plotted in FIG. [2l We 
associate the operators i?, Di and D2 to an empty lat- 
tice site, a first-class particle and a second-class particle 
respectively. These operators satisfy the algebra (jH]) for 
= 1,2. Their matrix representation is also given by ([3]) 
in which /^(i) = n}=i ITIJ) M = 2. It is much eas- 
ier to work in a grand-canonical ensemble in which the 
number of first-class particles is controlled by a fugacity. 
According to the matrix product method discussed above 
and following [11] the grand-canonical partition function 
of this system can be written as: 

Zn{z) = Tr{C^-^D2) (7) 

in which C — zDi + E. Now one can rewrite ([7]) using 
the matrix representation of D2 given in ^ as: 

N-l 

^ivW-E/2(j)(j|C^-'lO) (8) 

in which we have used the fact that: 

C\j)^zf,{j)\0) + \j + l). (9) 

Note that C^^^ on ]0) can be written as a linear super- 
position of {|0), • • • ,\N- 1)}. Expanding dH) as a power 
series in z, the canonical partition function of the gener- 
alized driven-diffusive system is equal to the coefficient 
of z^-\ 

This disordered driven-diffusive system can be mapped 
onto a HZRP consisting of M boxes and TV — AI sim- 
ilar particles. The boxes are divided into two different 
types. The rate at which a particle leaves the first box 
/i = 1 is given by ui(ni) while the rate at which a particle 
leaves the other boxes /i = 2, • • • , M is given by U2(n^). 
The mapping guaranties that the partition function of 
the disordered driven-diffusive system is equal to that of 
the HZRP. 

In [11] it has been shown that a much simpler HZRP in 
which the particles leave the boxes at a constant rate can 
be mapped onto a two-dimensional walk model. In that 
which follows we will show that the HZRP defined above 
can also be mapped onto a walk model. One can see 
that the expression ([5]) can be interpreted as the partition 
function of a simple random walk model if we consider C 
as the transfer matrix of a two-dimensional random walk 
defined on the lattice sites of a two-dimensional plane 
= ■ i,j > are integers}. Each step taken 

by the random walker has a weight. The random walker 
moves according to the following rules: 

• For i > j from the lattice site (i, j) to (i + 1, j + 1) 
with a weight 1 where i, j — 0,1,2, ■ ■ ■ , 00. This is 
called an upward step. 
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FIG. 2: A disordered driven-difFusive system on a ring of 
length A'' = 16 containing M = 5 particles. The first particle 
(white disk) is considered as a second-class particle while the 
rest of the particles (black disks) are considered as first-class 
particles. 

• For i > j from the lattice site to {i + 1,0) with 
a weight zfi(j) where z, j = 0, 1, 2, • ■ • , oo. This is 
called a horizontal (for j = 0) or a downward (for 
j ^ 0) step. 

Keeping in mind that the random walker does not take 
any steps in the negative z-direction, we assume that it 
starts from the origin (0, 0) and takes — 1 consecutive 
steps according to the above mentioned rules. Finally 
the random walker should be at the lattice site (A^ — 1, j) 
where j — 0, - ■ ■ ,N—1. We define the weight of a path of 
length A^ — 1 as a product of the weights of the steps taken 
by the random walker times an extra weight associated 
with the destination point that is /2(i) if the random 
walker's destination lattice site is (A^ — It is now 

easy to see why the transfer matrix of the walk model is 
exactly equal to C defined in ([7]) and dH). The partition 
function of the walk model is also given by ([8|). 

V. THE CANONICAL PARTITION FUNCTION 

The disordered driven-diffusive system defined in pre- 
vious section contains exactly M particles. When we map 



this system onto a HZRP it contains M boxes. As we 
saw the grand-canonical partition function of the disor- 
dered driven-diffusive system (or equivalently the HZRP) 
is equal to that of the walk model defined by the rules in 
previous section. It is easy to verify that if we are inter- 
ested in the canonical partition function of the disordered 
driven-diffusive system, we need the sum of the weights 
of the paths which contain exactly N — M upward steps 
or equivalently M — 1 downward and horizontal steps. 
This is because each downward or horizontal step con- 
tains one z in the weight of that step. In what follows we 
show how we can calculate the canonical partition func- 
tion of the disordered driven-diffusive system. 
The expression ([S]) can be rewritten as: 

Af-l 

Zn(z) = - ^0 - l)^(^o) (10) 

feo=0 

in which we have defined F(/co) = (OlC*^" |0). Clearly we 
have F{0) — 1. Using the matrix representation ([21) one 
finds that F(fco) satisfies the following recursion relation: 

feo-l 

F{ko) = J2 ^hih)F{kQ - fci - 1) (11) 

fci=0 

for /co = 1, 2, • • • , A^ — 1. After some calculations the 
recursion relation pip can be solved and one finds: 

^(fco)=E[n( E ^/i(%-i-%-i))]^/i(fc.-i)- 

j=0 j = l kj = l 

(12) 

Care should be taken in using in this compact form 
since its first term for z = is defined to be z/i(fco — 1). As 
we mentioned above, we are interested in the canonical 
partition function that is the coefficient of z^~'^ in (|10p . 
This can be achieved using ([TU)) and and the result 
is: 



Znm = E [ n ( E - - l))]/i(^M-2 ~ l)f2{N -ko- 1). (13) 

fco = l j = l fej = l 

I 



This is the canonical partition function of three differ- 
ent models which can be mapped onto each other; a 
disordered driven-diffusive system, a HZRP and a two- 
dimensional walk model as described in previous sections. 
Let us investigate two special cases studied in and 



[13]. In the first case we have Mi(7z) = 1 and U2{n) — p 
(equivalently /i(n) = 1 and /2(r^) — 4r). Using the fact 
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that: 



VI. CONCLUDING REMARKS 



EE- E 1 



fel=l fe2 = l fcjvf-2 = l 



fco-1 
M-2 



(14) 



in which 



b J ~ bi.ia-by. 
cient, it is easy to see that (IT5|) results in 



is the usual binomial coeffi- 



zn,m — 

ko=0 



N -ko-2 
M - 2 



(15) 



which is exactly the canonical partition function of the 
system studied in [l3| • In the second case introduced and 
studied in [l^ we have ui{n) = 1 and U2(n) = p(l + ^). 
For A = it converges to the first case mentioned before. 
It is easy to verify that fi{n) = 1 and /2(n) = p^j^^^^^ 
in which (a)„ is the Pochhammer symbol defined as: 

(a)„ =a(a + l)(a + 2)---(a + n-l) , (a)o = 1. (16) 

Replacing these into one finds: 



E 

fcn=0 



N -ko-2 
M -2 



(17) 



In terms of the number of empty lattice sites, that is 
L := N — AI, one recovers the expression Zl,m which 
has already been obtained in [l7j . 



In this paper we have investigated a disordered driven- 
diffusive system defined on a ring which can be mapped 
onto a HZRP. We have shown that the grand-canonical 
partition function of the HZRP can be interpreted as 
the partition function of a two-dimensional walk model. 
We have obtained exact expressions for both the grand- 
canonical and canonical partition functions of the HZRP. 
We have presented two simple examples which confirm 
the correctness of our calculations. Apart from the 
special driven-diffusive system studied in this paper, in 
which a single second-class particle beside a finite number 
of first-class particles, lie on the lattice and hop stochasti- 
cally in a preferred direction, one can also consider other 
cases. For instance, a disordered driven-diffusive system 
in which the backward hopping of the particles is also in- 
cluded similar to the one studied in [lj| . In terms of the 
walk model this results in horizontal movements at all 
heights. In a variety of the zero-rage processes have 
been investigated. It would be interesting to investigate 
the walk models associated with different types of these 
processes. 
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